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Abstract 

We study the phase behavior of diblock copolymers in presence of an external electric field. We 
employ self-consistent field theory and treat the relevant Maxwell equation as an additional self- 
consistent equation. Because we do not treat the electric field perturbatively, we can examine its 
effects even when its magnitude is large. The electric field couples to the system's morphology only 
through the difference between the dielectric constants of the two blocks. We find that an external 
field aligns a body-centered cubic phase along the (111) direction, reducing its symmetry group 
to R3m. Transitions between this phase and the disordered or hexagonal phases can occur for 
external electric fields ranging from a minimum to a maximum value beyond which the R2m phase 
disappears completely. This electric-field range depends on diblock architecture and temperature. 
We present several cuts through the phase diagram in the space of temperature, architecture and 
applied field, including one applicable to a system recently studied. 
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I. INTRODUCTION 



Because block copolymers readily self-assemble into various ordered arrays, they have 
been avidly studied for technological applications such as high density porous materials, 
nano-lithographic templates, photonic band gap materials 

HQS 

, and well-ordered arrays 

of metal nano- wires 1, ^J. One practical difficulty to their use in some applications is that 
the ordered phase is not created in one single crystal, but rather in domains of differing 
orientation. One means of aligning the domains is to apply an external electric field. It 
has been shown that applying an electric potential, on the order of a few to 

a few dozen volts across electrodes separated by several micrometers, can effectively orient 
domains of lamellar or cylindrical morphology normal to the surfaces of thin films. The basis 
of this orientation effect is simple. In order to reduce accumulation of polarization charge, 
the system lowers its free energy by aligning cylinders or lamellae so that their long axis is 
parallel to the applied field. 

Recently, related experiments on diblock copolymers JjJ have been performed where ex- 
ternal fields have been applied to bring about a phase transition from a phase of spheres to 
one of cylinders. In the phase of spheres, it is not possible to eliminate the accumulation 
of polarization charge so that its free energy increases in an external field with respect to a 
cylindrical phase, and a phase transition can be induced. This change in phase due to the 
application of an electric field was considered by Tsori et al. [UJ, and by Xu et al [9]. 

The effects of an external field on an ordered array of inhomogeneous dielectric material is 
of great interest. First, the problem is inherently self-consistent simply because the material 
is a dielectric; i.e. the electric field at a given point depends upon the polarization at that 
point which, in turn, depends upon the local electric field. In addition, in the problem of in- 
terest here, the local dielectric constant is inhomogeneous. It depends upon the morphology 



111- 



of the ordered phase, which itself depends upon the local electric field 

In previous calculations for diblock copolymers 0,0], this self-consistent circle has been 
broken by assuming that the two block are only weakly segregated, resulting in a small 
amplitude of the spatial variation of the relative concentration of the two blocks. In this case 
it follows from the vanishing of the divergence of the electric displacement that the amplitude 
of the spatially varying electric field is also small so that the electrostatic Maxwell equation 
can be solved perturbatively for the electric field as a function of the order parameter. 
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This procedure was carried out to quadratic order in the field by Amundson et al.|l2j. It 
is appropriate in the weak-segregation limit, and should be adequate for determining the 
general phase behavior in weak external fields. However, since experiments are often not in 
the weak segregation limit, and the effect of electric fields has hardly been explored, further 
study is clearly called for. 

In recent years, thermodynamic properties of block copolymer systems have been treated 
successfully by the full self-consistent field (SCF) theory, to which weak- and strong- 
segregation theories are approximations [13]. Given the self-consistent nature of an in- 
homogeneous dielectric in an external electric field, it seems natural to apply the full SCF 
theory to this problem as well. That is what we do in this paper. We solve exactly the 
full set of SCF equations and the appropriate Maxwell equations under the assumption of 
a simple constitutive relation between the local dielectric properties and the local volume 
fractions. In particular, we consider the evolution of the bulk phase diagram of diblock 
copolymers in an applied electric field, and focus upon its effect on reducing the region 
of the phase diagram occupied by the body-centered cubic (bcc) structure, (space group 
7m3m). Evolution of the gyroid structure (space group Ia3d), whose region in the phase 
diagram also decreases due to the application of a field, is not considered. 

We calculate the strength of an external field needed to bring about a phase transition 
from the (distorted) spherical phase to the disordered phase and to the cylindrical phase. 
For the transition to the latter phase we find two distinct behaviors depending upon the 
architecture of the diblock, as measured by the parameter fa introduced below. The first is 
brought about if a transition from the spherical phase to the cylindrical phase can be induced 
in the absence of an external electric field simply by reducing the temperature in the realm 
of interest. If so, the same must also be true for very small fields. As a consequence, one 
can always find a temperature in that realm at which an arbitrarily small field will induce 
a transition from the spherical to the cylindrical phase. The other behavior occurs if the 
spherical phase is the most stable one in the absence of an external field for temperatures in 
the realm of interest. In that case, a non-zero external field is required to induce a transition 
from it to the cylindrical one at any temperature in this realm. In either case, we find that 
for a given architecture, there is a maximum value of applied field beyond which the spherical 
phase is no longer the most stable one for any temperature. 

In the following section, we set up the general formalism. In section III, we discuss 
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its application to the phase of (distorted) spheres, and compare the results of the full self- 
consistent calculation with those obtained from an expansion of the free energy in the electric 
field to order E 2 . Such an expansion does not indicate the optimal direction in which the 
field aligns the cubic phase, whereas the full calculation shows that alignment along the 
(111) direction is favored over a (100) orientation. There is a concomitant reduction of the 
symmetry of the phase from Im3m (bcc phase) to R3m (distorted spherical phase). Various 
cuts through the phase diagram are also presented. We conclude with a brief summary and 
comparison with recent experiments. 



II. GENERAL FORMALISM 



We consider a melt of n A-B diblock copolymer chains, each of polymerization index 
N = Na + Nb- If the specific volumes of the A and B monomers are va and vb, respectively, 
the volume per chain is v p = NaVa + NbVb- For an incompressible melt of A-B chains, the 
volume fraction of the A monomers is NaVa/ {NaVa + NbVb), and the total system volume 
is Q = nv p . We assume the monomer volumes to be identical, va = vb, so that the volume 
fraction of the A-monomers is equal to the mole fraction of the A-monomers, Ja = Na/N. 
We also assume that the Kuhn lengths of the A and B components are identical, a length 
denoted a. 

In the absence of an external field, the application of SCF theory Q| leads to a free energy 
T which is a functional of unknown fields Wa, Wb, and S, and a function of temperature T 

T(Wa Wr E-T) 1 r 

1 . ' ' ^ = - In Q[W A , Wb} + jz dr { X N$ A $ B ~ W a $a ~ W b $b 
nk B T \l J 

-S(l-$ A -$ fl )} , (1) 

where ks is the Boltzmann constant, &a{y) and $_e(r) are the local volume fractions of A 
and B monomers. The dependence on T comes from the usual Flory interaction parameter, 
X, which to a good approximation is inversely proportional to the temperature, xN = b/T 
with b a constant. The function Q[Wyi, Wb] is the partition function of a single polymer 
chain subject to the fields W^(r) and Wb{t), as is given below. The field S(r) is a La- 
grange multiplier that enforces locally the incompressibility constraint, ^a{ y ) + $s( r ) = !■ 
The three unknown fields are determined by requiring that the free-energy functional be 
extremized with respect to their variation at constant T. 
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The fields W A and Wb appear in the single-chain partition function of the flexible diblock 
copolymer, Q[Wa, Wb] = J drq(r, l)/c, where q(r, s) satisfies the modified diffusion equation 

|j = l -Na 2 V 2 q - W A {v)q, if < s < f A , (2) 



and 

0(1 l Na 2 V 2 q-W B {v)q, if f A < s < 1 , (3) 



ds 6 

with the initial condition q(r, 0) = 1, and c is a volume of no consequence here. 

The addition of a local electric field E(r) in the derivation of the free energy T is straight- 
forward. In an ensemble for which an external electric potential is held fixed Q], the above 
free energy simply becomes 

+ n I dr {* mA ® B ~ Wa ® a ~ Wb ® b ~ s(1 -®a-§b)} , (4) 

where eo is the vacuum permittivity, and n(r) is the local dielectric constant. A constitutive 
relation between n(r) and the volume fractions of A and B monomers must be specified. We 
choose a linear interpolation relation 

k(t) = K A $4(r) + K B <£> B (r) , (5) 

where k a and Kb are the dielectric constants of pure A and B homopolymer phases, respec- 
tively. This choice is clearly correct in the limiting cases of the pure systems, and in the 
weak-segregation limit. It also has the virtue of simplicity and should capture the correct 
physics. 

From the above it can be seen that a convenient scale for the strength of the electric field 

is 

£ =(—) ■ (6) 
The magnitude of this electric field unit at typical experimental temperatures, T ~ 430 K, 
and for typical volume per polymer chain, v p ~ 100 nm 3 , is £ ~ 82 V//zm. We shall denote 
the dimensionless electric field rescaled in this unit as E = E/£. Similarly a dimensionless 
displacement field, D, is conveniently defined by D = D/eo<£. 

The requirement that the free energy functional be an extremum with respect to variation 
of W A , Wb, H, and of the volume fractions $> A and $b at constant temperature, or %iV, and 
fixed electric field E, leads to the following set of SCF equations: 



w A = X N<p B + !i- ^a|E| 2 , (7) 
w B = X N0 a + £-^k b \E\ 2 , (8) 
(Pa + <t> B = 1 , (9) 

*■ " "If • 

The values of Wa, Wb, H, $a and $g, which satisfy these equations are denoted by lower 
case letters, wa, Wb, £, <Pa, <Pb, respectively. The free energy within the SCF approximation, 
F sc f is obtained by substitution of these values into the free energy of Eq. (p|J), 

F scf (T, E) = F(w A , w B , £; T, E) , (12) 

or 

= - In Q[w A , wb]~^J dr{ X NMr)Mr) + £(r)] • (13) 

In addition to these equations, there are also the Maxwell equations which the electro- 
static field must satisfy in absence of free charges: 

V x E = , (14) 
V • D(r) = V • (e„K(r)E(r)) = . (15) 

As usual, we guarantee that the first of these equations is satisfied by introducing the electric 
potential V^(r), 

E(r) = -W(r) = -W(r)/e . (16) 

Since we will consider, in addition to the disordered phase, spatially-periodic ones, it is 
convenient to write all functions of position in terms of their values averaged over a unit cell 

C = (C) = •^ccuC'(r) dv ^ (17) 

Junit cell <^ r 

and their deviations from those average values 

5C(v) = C(r) - Co . (18) 
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The average values of several quantities of interest are 



4>A,0 = fA , 
4>B,0 = 1 - fA , 

«0 = K A f A +K B (1- J'a) , 

wa,o = XN(1 - f A ) - ^|E | 2 , 

w B ,o = xNfA-lK B \E \ 2 , (19) 

where the value of £o nas been arbitrarily set to zero, and E is the value of the local 
electric field averaged over a unit cell. To determine this without knowing the full spatially 
dependent electric field E(r), we reason as follows. Assume that the external field is produced 
by planar electrodes which are separated by a distance d and subject to a voltage difference 
Vi2- In the gap, and along the z-axis perpendicular to the electrodes, the field is E ext = 
—V12/L. Given that the dielectric fills the space between the plates, and that the voltage 
V12 is held fixed as the dielectric is inserted, it follows that J E z dz = E ext L, and that the 
average value of E z is E ext . We make a reasonable assumption that the free energy of the 
system is minimized when an axis of symmetry of one of the ordered structures coincides 
with the z-axis. In this case Eq = J L E z dz/L = E cxt . Hence in rescaled units 



En (^V /2 £cxtZ, (20 ) 



k B T 
and 

5E(r) = E(r) - E = -V<5V(r) . (21) 

Utilizing these average values, we can rewrite the free energy in the SCF approximation, 
Eqs. (TT2j) -(fB J) . in the form 



-T% = " ' n ( o? K ' mB M + ~ W " \-°El - *g- I S M r)SMr)dr , (22) 

nk B T {Q{w A fl,w Bfi \} 2 il J 

where, from the incompressibility condition, 54>a{t) = —8<ft B (r). Note that the electric- field 
contribution —k Eq/2 is common to all phases. For the lamellar and hexagonal phases in the 
lowest energy orientation, this is the only contribution to the free energy from the electric 
field. It can conveniently be absorbed in a redefinition of the free energy 

f n (E Q ) = E sci /nk B T + ^k E 2 . (23) 



The advantage of separating out the average values is that the only remaining Maxwell 
equation, Eq. (|15p. can be written as an inhomogeneous equation for the potential SV(r), 

V8V(r) ■ V[K A 8<J) A (r) + kb50b(f)] + [«a(/a + S<j> A (r)) + «b(1 - /a + ^(r))] V 2 5V(r) 
= £ ^A(r) + K B 5(j) B (r)] • (24) 
This, with the three remaining self-consistent equations, 

5w A (r) = X N6M*) + ^(r) + ~k a [2E • V6V(r) - (V5V(r)) 2 ] , (25) 

5w B (r) = xN6<j> A (r) + <^( r ) + 7^b[2E • V8V(r) - (V5V(r)) 2 } , (26) 
5<p A {v) + 50 B (r) = , (27) 

constitute the four self-consistent equations which determine the four functions 5w A (r), 
5w B (r), 8Z(r), and 8V(r). 

We note that with our choice of constant external field applied along the z direction, the 
Maxwell equation, Eq. ()24jl . admits the following symmetry: 



SV(r L ,z) = SV(-t±,z) = -SV(t ± ,-z) 
n(r ± ,z) = n A (f) A (r ± ,z) + K B (p B (r ± ,z) 

= «(-rjL,z) = k(t±,-z) , (28) 

where the components of r have been written as (r±,z). The self-consistent equations, 
Eq. (j2I| -()27jl. are now solved by a standard procedure of expanding the functions of position 
in a complete set of functions with the above symmetries and those of any specific phase 
considered We have utilized in our calculation sets of basis functions containing between 
70 and 125 functions, depending upon the value of %iV '. 

The only parameters entering our calculation are \N ■, /a, &a and k b , and the rescaled 
external field E cxt /£- Comparison of the results with experiment requires the evaluation of 
8 for given T and volume per chain v p . In addition, the relation between T and xiV must 
be specified. 



III. RESULTS 



As noted earlier, the free energies of lamellar or hexagonal phases are minimized when the 
lamellae or the cylinders are aligned parallel to the electric field because in this orientation 
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there is no buildup of polarization charge. In the body-centered cubic (bcc) phase, however, 
there must be an accumulation of polarization charge irrespective of the external field direc- 
tion. We must determine which field direction produces a phase of distorted spheres with 
the lowest free energy, and how large a field this phase can sustain before a transition to the 
hexagonal phase is encountered. It is these issues which we now address. 



A. The R3m to Hexagonal phase transition 

The symmetry group of the bcc phase is Im3m which has three two-dimensional space 
subgroups: p4mm along the [100] direction, p6mm along the [111] direction, andp2mm along 
the [110] direction. If the field were applied along either the [110] or [100] directions the 
symmetry would be reduced to I4/mmm, while if it were applied along the [111] direction, 
the symmetry would be reduced to R3m. The symmetry in the latter case is of a bcc 
arrangement of spheres that has been distorted along the [111] direction. As the R3m group 
has the p6mm symmetry of the hexagonal phase, one would suspect that a field applied 
along the diagonal [111] direction will result in the lowest free energy. By direct calculation 
of these configurations, we find that the R3m phase does indeed have a lower free energy 
than that of the IA/mmm. 

That the electric field favors one orientation of the Im3m over another is an effect which 
is not captured by an expansion of the free energy to quadratic order in the external field 



111 ] . Nonetheless it is instructive to consider the result of such an expansion. It is obtained 
by solving the Maxwell equation V ■ [e K(<p A , 0s)E] = to second order in E to obtain 
E(4>a, <Pb), and evaluating this field from the volume fractions characterizing the system in 
the absence of an external field. The distortion of the density distribution produced by this 
field itself contributes terms to the free energy which are higher order in E 2 . For a phase 
which is cubic in the absence of an electric field, the perturbation result can be written as 
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F pt (Eo) 
nk B T 



--k cS Eq, 



1 - 



1 f K A 



12Q. V 



Kq 



dr[S(f> A (r) - 5<p B {v)f 



(29) 



where 



is the variation of the local volume fraction in the zero-field structure, 
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and K e g is, by definition, the effective dielectric constant for the structure in the field. 

We now compare the full SCF solution with this perturbation result. We choose \N = 15 
and Ja = 0.29, values at which the bcc phase is the most stable in zero electric field. 
The dielectric constants are chosen to make contact with recent experimental systems of 
polymethylmethacrylate (PMMA) /polystyrene (PS) diblock copolymer, which is referred 
to hereafter as the PMMA-PS system. At experimental temperatures around 160° C the 
dielectric constants appropriate to the PMMA-PS copolymer with PMMA being the A 
block and PS the B block are: k a = 6.0 (for PMMA), and k b = 2.5 (for PS) which 
yield an average of k,q ~ 3.52. In Fig.^ we show the difference, Af n , between the free energy 
f n (Eo) = F sc f/nkBT + k,qEq/2, Eq. (|23|) . and its value in zero external field in the bcc phase. 
It is shown as a function of Eq, for the hexagonal phase and for the R3m phase, as calculated 
from the full SCF theory and from perturbation theory. The latter is seen to be adequate 
for fields smaller than ten to twenty percent of the natural unit £ at which E = 1. The 
figure also shows that there is a transition from the R3m to the hexagonal phase at a value 
of Eq ~ 0.477 as determined from the full self-consistent calculation. Perturbation theory 
underestimates the magnitude of the field needed to bring about this transition. That the 
transition is first-order is easily seen as follows. The average electric and displacement fields, 
Eq and Do, are evaluated by taking their spatial averages over the unit cell. In our case the 
only non-zero average components are those in the z-direction, and they are related to the 
free energy per unit volume according to 

dF/Q 



OEq 

or 



■D , (30) 



^ = ■ (3D 
Ohio 

One sees from Fig.^that at the phase transition, the free energies of the R3m and hexagonal 
phases intersect with different slopes, therefore the displacement field changes abruptly. 

As a result of the application of the electric field along the [111] direction, the spheres 
of minority component are elongated in this direction. A density profile of the system in 
the R3m phase at an external field Eq = 0.470, slightly smaller than that at the transition 
to the hexagonal phase, Eq = 0.477, is shown in Fig. Efb). At the transition, the profile 
changes abruptly to that of the hexagonal phase, which is also shown in|2fc) for Eq = 0.480. 
To see the extent of the distortion in the R3m phase, which can be characterized by the 
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aspect ratio of the distorted spheres, 1.248, we also present the density profile of the bcc 
phase in zero external field, in|2(a). The cuts are in the plane containing the [111] and [110] 
directions. 

There are two features of interest that can be seen particularly clearly from the approx- 
imate expression of Eq. (|29|) . The first is that in the R3m phase, the effective average 
dielectric constant, K e g, is smaller than kq. In the hexagonal and disordered phases, how- 
ever, K e g is precisely Kq- Therefore, the displacement field Dq in the R3m phase is smaller 
than in the other two phases. This is in accord with the change of slope of the free energy 
with electric field shown in Fig. [T] and Eqs. (|BT)J) - (fB*T]) . 

The second concerns the fact that the dielectric constants are temperature dependent. 
Therefore, the value of the electric field needed to bring about a phase transition will also 
vary with temperature. The perturbation expression leads one to expect that, for fields 
smaller or comparable to S, the natural .E-field scale, the field at the transition will vary as 

p m „ MO] 1/2 [f A KA(T) + (l-f A )K B (T)]V* 

tr{ } « A (T) - k b (T) ~ k a {T) - k b {T) ■ {62) 

B. The generalized Claussius-Clapeyron equation 

Before presenting the phase diagram of our A/B block copolymer system in an E-Held, 
we will make use of some general thermodynamic considerations. In particular, from the 
differential of the free energy per unit volume 

d(F/n) = -sdT - D dE , (33) 

where s = S/Q is the entropy per unit volume, one immediately derives a Claussius- 
Clapeyron equation for the slope of the coexistence line between any two phases 

^ = (34) 
dT AD ' 1 ' 

where As and AD are the differences in entropies and displacement fields, respectively, in 
the coexisting phases. This can be expressed in terms of Eq, Dq and xN = b/T as 

dE v p A(s/k B ) E 



d( X N) X N AD 2 X N ■ 
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C. Phase diagrams 

We now turn to the phase diagram as a function of the inverse temperature, xiV, the 
A-monomer fraction, Ja, and the applied external field, E = E /£. We concentrate on the 
portion of the phase diagram involving the phase with R3m symmetry and the neighboring 
disordered and hexagonal phases. In the space of inverse temperature ^iV, the fraction 
/a, and applied field, the R3m phase occupies a volume which is bounded by two sheets 
of first-order transitions: one from the R3m to the hexagonal phase, the other from the 
R3m to the disordered phase. These two sheets of first-order transitions meet at a line of 
triple points, [-^.triple (/a) > X-^tripie(/A)]- Beyond this line, the R3m phase no longer exists, 
while the disordered and hexagonal phases remain. They are separated by another sheet of 
first-order transitions which emerges from the line of triple points. Hence this line is the 
locus at which all three sheets of first-order transitions meet. 

In Fig. |3] we show a cut through the phase diagram at fixed A-monomer fraction, Ja = 
0.29. The cut shows the phase diagram as a function of the dimensionless electric field Eq and 
xN. At zero external field, the entropy difference between the bcc phase and the hexagonal 
phase is non-zero, but the difference in displacement field obviously vanishes. From the 
Claussius-Clapeyron equation, Eq. (J33J), the slope of the phase boundary between these two 
phases must be infinite at zero E-field. The same is true for the slope of the phase boundary 
between the bcc phase and the disordered phase at vanishing E-fields. Furthermore, we 
know from the zero electric field results that the entropy of the disordered phase is greater 
than that of the bcc phase which, in turn, is greater than that of the hexagonal phase. We 
also know that the displacement field in the disordered and in the hexagonal phases is equal 
to KoeoEo. As we noted earlier, the displacement field in the R3m phase is less than this 
value. This information, together with the Claussius-Clapeyron, Eq. (|35jh implies that the 
phase boundary between R3m and the disordered phase has a negative slope, while that 
between R3m and the hexagonal phase is positive in accord with Fig. El 

Moreover, because of the presence of the positive second term in Eq. (J35)) . the positive 
slope of the phase boundary between disordered and R3m phases will be greater in mag- 
nitude, or steeper, than that between the R3m and hexagonal phases. This is borne out 
by Fig. El The three phases meet at the triple point, above which the phase boundary is 
vertical as there is no difference between the displacement fields of the coexisting disordered 
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and hexagonal phases. The value of the electric field at the triple point is £o,tripie — 0.71. 

To make contact with experiment, we take parameters to fit the PMMA-PS system of Ref 
Pi . With j a = 0.29 and a molecular mass of 3.9 x 10 4 g/mol, and utilizing the known values 
of monomeric volumes, we obtain a chain length of iV ~ 379 and a volume per PMMA-PS 
chain of v p = 61.24 nm 3 . At T = 430 K this yields 8 = (k B T/e v p )^ 2 = 104.6 V///m. 
Therefore, the value of the electric field at the triple point is in physical units £o,tri P ie ~ 74.5 
V//xm at this value of Ja and T. One sees from the figure that a transition from R3m 
to hexagonal phases could be brought about at electric fields within the interval from this 
maximum value down to zero, depending upon the values of \N and Ja- 

The evolution of the phase diagram of Fig. E] with A-monomer fraction, Ja, is easily 
understood. As Ja decreases from 0.29, the phase boundary at zero field between R3m 
and hexagonal phases moves toward greater values of xN as does the boundary between 
disordered and R3m phases. When Ja is smaller than /^ ocx = 0.114, the value at which the 
bcc and hexagonal phases coexist at infinite xN [ IS] , the boundary between hexagonal and 
R3m phases will asymptote with zero slope to an /^-dependent finite value as xN increases 
without limit. This zero slope also follows from the Claussius-Clapeyron equation ()35|) due 
to the fact that the ratio As/AD is finite and 1/xiV — > 0. 

An example of such a phase diagram is shown in Fig. 0] This figure corresponds to a 
system with J'a = 0.1 < /™ ex = 0.114, as was investigated recently in Ref. In contrast 
with Fig. El one sees here that the interval over which a transition can be observed from 
R3m to hexagonal phases now extends from the triple point at Eq = 2.56 down to a non-zero 
minimum value of E = 1.33. That is, for electric fields less than this minimum value, no 
transition from the R3m to a hexagonal phase occurs within our model. For PMMA-PS with 
/a = 0.1, and molecular mass of 1.51 x 10 5 g/mol, as in Ref. |9(, we obtain N rs 1458 and a 
volume per chain v p = 239. 7 nm 3 . Therefore at the experimental temperature of T = 430K, 
the unit of electric field 8 = 52.9 V//im. In physical units, then, the triple point occurs at 
an external field of about 135V//xm and the minimum electric field needed to produce a 
transition can be estimated to be 79 V//zm. 

In Fig.EJwe show a different cut through the phase diagram in the (f?o> Ja) plane and for 
a fixed xN — 13.3. The location of the triple point is -Eo,tri P ie = 0.58 and /A,tri P ie = 0.320. 
This figure, and that of Figs. El and 01 show that the value of the electric field needed to 
bring about a transition from the R3m phase is, for a given Ja fraction, a sensitive function 
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of temperature and, for a given temperature, a sensitive function of the mole fraction of A 
block, f A . 

Figure shows two cuts through the phase diagram at constant electric field, E = 0, 
and Eq = 0.2. The solid line at lower xN shows the phase boundary at zero-field between 
disordered and bcc phases while the solid line at larger \N shows the zero-field phase 
boundary, between the bcc and hexagonal phases. The dashed lines between them show the 
phase boundaries for Eq = 0.2. The line denoted B is the boundary between the disordered 
and the R3m phase of distorted spheres, A is the boundary between R3m and hexagonal 
phases. These boundaries meet at the triple point, tr, which occurs at xiV « 11.43, and 
] a ~ 0.39. For larger values of f A , there is a line, C, of transitions directly from the disordered 
to the hexagonal phase. As the external field increases still further, the triple point recedes 
to larger values of \N leaving behind only the line of direct transitions between disordered 
and hexagonal phases. Note that, except for the location of its terminus at the triple point, 
this boundary is independent of the applied field as it contributes to the free energy of both 
of these phases equally. The dielectric constants used to generate this figure are the same 
as those used in previous figures. 

For completeness, we have also examined the case in which the dielectric constants of 
the minority and majority components are interchanged as compared with Fig. El Namely, 
the majority component with, f A = 0.71 has the larger dielectric constant of ka = 6.0 and 
the minority the smaller value of Kg = 2.5. We find that the R3m phase is now somewhat 
more stable with respect to the hexagonal phase, so that the value of the external electric 
field needed to bring about a transition from the former to the latter phase is increased. We 
note that this interchange increases the average value of the dielectric constant, so that all 
phases have a lower free energy due to the factor of —KqEq/2 which it contains. However, 
it is not a priori obvious that the R3m phase would have its free energy lowered by more 
than that of the hexagonal phase by this interchange. In addition, we have determined 
that the spheres of minority component and lower dielectric constant distort in the [111] 
direction just as in the case when the minority component has the larger dielectric constant. 
The above effects are not captured by the perturbation result of Eq. (}2*9"|) which is invariant 
under the interchange of ka and Kg. 
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IV. CONCLUDING REMARKS 



In sum, we have calculated the phase diagram of a block copolymer system in an ex- 
ternal electric field which couples to the diblocks through the difference in their dielectric 
constants. We have employed a fully self-consistent field approach in which the relevant 
Maxwell equation is treated on an equal footing with the other self-consistent equations. 
We have determined that the body-centered cubic phase will preferentially align along the 
[111] direction causing its symmetry to be reduced to R3m. The electric field can induce 
phase transitions between this phase and either the disordered or the hexagonal phase. The 
strength of the field needed to induce such transitions is a sensitive function of the parame- 
ters of the system, such as its temperature and its chain architecture, which in the case of 
linear diblocks is quantified simply by the mole fraction, J'a- 

For parameters that fit the experimental PMMA-PS diblock copolymer system investi- 



gated recently jjj: Ja = 0.1, v p = 239.7nm a , and T = 430 K, we find that an electric field of 
at least 70 — 80 V//xm would be needed to observe a transition to the hexagonal phase. This 
contrasts with the reported existence of such a phase transition under an applied field of only 
40 V/ /iin. There are several possible explanations of the difference between the experimental 
results and the theoretical ones presented here. 

Our model employs a linear constitutive relation between dielectric constant and volume 
fractions, and characterizes the PMMA-PS system by a few general parameters, the PMMA 
mole fraction Ja, and the interaction parameter It further assumes equal volumes for 
both monomers and equal Kuhn lengths for them. One knows that deviations from the 



last assumption certainly shift the locations of the phase boundaries [17 1. It is plausible 
that at rather asymmetric volume fractions of J'a = 0.1, the model provides only semi- 
quantitative agreement with the experimental PMMA-PS phase diagram. Any difference in 
the theoretical and experimental phase diagrams at zero electric field will, in the presence of 
a non-zero one, manifest itself in a difference in relative stability of the various phases. Given 
the sensitive dependence on the phase boundaries of the minimum external field needed to 
bring about a phase transition, differences between the general theory and the experimental 
result are to be expected. At present, the phase diagram of the PMMA-PS system of Ref. 
is not yet known. When additional experimental information becomes available, one will 
also need to determine the relationship between the temperature and the xiV interaction 
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parameter in order to convert the phase diagram calculated here to practical units so that 
it can be compared directly to the experimental one. 

Lastly, we have employed a simple coupling between the system and the external field via 
the difference in dielectric constants of the copolymer blocks. Other couplings are possible 
[isL Jioil • Just such an additional coupling, to mobile ions, has been suggested by Tsori et 
al. [10| and is discussed also in Ref. |9(. A minute fraction of mobile ions embedded in the 
minority PMMA fraction and not in the majority PS can lead to an enhanced response of 
the PMMA-PS system to external electric fields with moderate magnitude. It could also 
change the phase diagram quantitatively resulting in a substantial lowering of the triple- 
point value of the electric field. Additional experiments, particularly on copolymers with 
the same PMMA-PS blocks, but at different temperatures or values of the architectural 
parameter Ja, would be most useful to shed additional light on the comparison of theory 
and experiment. In particular, a comparison of the two PMMA-PS systems of Ref. jj| and 
Ref. p would be enlightening because, as Figs. Eland lUshow, they are predicted here to 
exhibit significantly different phase behavior in an external electric field. 
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FIG. 1: The difference, A/ n , between the dimensionless free energy, f n , defined in Eq. (|23|) and its 
value in zero external field in the bcc phase. It is calculated from the SCF theory, and is plotted 
versus dimensionless external electric field, Eq = Eq/£, for the hexagonal phase (horizontal dotted 
line) and the R3m phase (solid line). The system is characterized by a = 15, f a = 0.29. The 
dielectric constants are: ka = 6.0 (for the PMMA block), and kb = 2.5 (for the PS block), yielding 
kq ~ 3.52. The perturbation theory result for the R3m phase is shown as a dashed and dotted 
line. It has a higher free energy. 
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FIG. 2: Density profiles for three different phases of a system characterized by \N = 15 and 
/ = 0.29, with other parameters as in Fig. 1. (a) the bcc phase which occurs in zero external 
field; (b) the R3m phase at an external electric field Eq = 0.470 just below the phase transition 
to the hexagonal phase which occurs at Eq = 0.477. (c) the hexagonal phase, which is shown 
for Eq = 0.480. The cuts are in the plane containing the [111] and [110] directions. In the black 
regions, the local volume fractions of component A is greater than 0.55, in the intermediate regions, 
it is between 0.55 and 0.45, and in the white regions, it is less than 0.45. 
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FIG. 3: Calculated phase diagram of a diblock copolymer with a volume fraction of /a = 0.29 
in the presence of an external electric field. The phase diagram is shown as a function of the 
dimensionless field Eq and the interaction parameter \N ■ The triple point is located at .Eb.tripie = 
0.71, %iV = 14.58 tr ipie- Other parameters as in Fig. 1. 
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FIG. 4: Calculated phase diagram of a diblock copolymer in the presence of an external electric 
field. Similar to Fig. [3 but with fraction of the A block, /a = 0.1. The phase diagram is shown 
as a function of the dimensionless field Eq and the interaction parameter %iV. The triple point is 
located at -Eo.tripie = 2.56, x^tripie = 50.74. 
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FIG. 5: Calculated phase diagram for a diblock copolymer as a function of dimensionless external 
field and the A mole fraction parameter, /a- Other parameters are \N = 13.3, ka = 6.0 and 
kb = 2.5. The triple point occurs at i?o,tripie = 0-58 and /U,tripie = 0.32. 
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FIG. 6: Calculated phase diagram at constant electric field, E . The outer two solid lines are 
the Eq = disorder-to-bcc and bcc-to-hexagonal phase boundaries. Between them we show three 
other transition lines for E$ = 0.2. They are the i?3m-to-Hexagonal, (A), i?3m-to-Disorder (B), 
and Disorder-to-Hexagonal transition (C). These three lines meet at tr, the triple point: /a. triple = 
0.390 and x^tripic = 11.43. 
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